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ABSTRACT 

The question of a purely hydrodynamic origin of turbulence in accretion disks 
is reexamined, on the basis of a large body of experimental and numerical evi- 
dence on various subcritical (i.e., linearly stable) hydrodynamic flows. 

One of the main points of this paper is that the length scale and velocity 
fluctuation amplitude which are characteristic of turbulent transport in these 

— 1/2 

flows scale like Rem , where Rem is the minimal Reynolds number for the onset 
of fully developed turbulence. From this scaling, a simple explanation of the 
dependence of Rem with relative gap width in subcritical Couette- Taylor flows 
is developed. It is also argued that flows in the shearing sheet limit should be 
turbulent, and that the lack of turbulence in all such simulations performed to 
date is most likely due to a lack of resolution, as a consequence of the effect of 
the Coriolis force on the large scale fluctuations of turbulent flows. 

These results imply that accretion flows should be turbulent through hy- 
drodynamic processes. If this is the case, the Shakura-Sunyaev a parameter 
is constrained to lie in the range 10~^ — 10~^ in accretion disks, depending on 
unknown features of the mechanism which sustains turbulence. Whether the 
hydrodynamic source of turbulence is more efficient than the MHD one where 
present is an open question. 

Subject headings: hydrodynamics - turbulence - accretion disks 



1. Introduction 



The need for turbulent transport to account for the rather short accretion/ejection 
time-scales of Young Stellar Objects (YSOs) and binary systems (CV, X-ray binaries), or 
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for the very large energy output of Active Galactic Nuclei (AGNs), is a well-known feature 
of accretion disk theory. From the very beginning, differential rotation has been regarded as 
one of the most promising sources for turbulence, since shear flows are known to be able to 
feed both hydrodynamic and MHD instabilities. 

In most accretion disk models, the angular velocity profile satisfies Rayleigh's criterion, 
implying that the corresponding hydrodynamic flow is linearly stable. This is the case 
in particular of the nearly keplerian velocity profile of cold disk models. However, finite 
amplitude instabilities are theoretically known to occur in some linearly stable fiows, and 
are believed to cause the turbulence observed in actual experiments, e.g. in planar Couctte 
fiows, or Couette- Taylor fiows with the inner cylinder at rest. Furthermore, shear- driven 
hydrodynamic turbulence would certainly produce the required outward transport of angular 
momentum for keplerian fiows, due to their outwardly decreasing angular velocity profile. 
For these reasons, turbulence in accretion disks (magnetized or not) has long been widely 
believed to originate in purely hydrodynamic phenomena. 

This picture has seriously been challenged in the past decade. First, Balbus and Hawley 
(1991) have shown that a local version of the magneto-rotational instability (Chandrasekhar 
1960) operates in differentially rotating disk. This instability was later recognized to give 

rise to MHD turbulence and transport as well as to magnetic field amplification (Hawley 
et al. 1995; Brandenburg et al. 1995). The physics of the magneto-rotational instability is 
by now a well-established aspect of accretion disk theory. Secondly, recent simulations of 
hydrodynamic fiuid fiows in the shearing sheet approximation strongly suggest that accretion 
disk fiows cannot become turbulent through hydrodynamic processes alone (Balbus et al. 
1996; Hawley et al. 1999); indeed, in these simulations, planar Couette fiows are observed 
to be turbulent, but turbulence disappears as soon as a Coriolis force is added, suggesting 
that this force prevents the onset of the finite amplitude instabilities through which linearly 
stable fiows are believed to become turbulent. Although this last finding seems to confiict 
with the available experimental evidence on Couette- Taylor fiows (Richard and Zahn 1999), 
it has strengthened the idea that linear magnetic instabilities play a key role in the onset of 
turbulence in accretion disks. 

The main objective of this paper is to critically reinvestigate the possibility of hydro- 
dynamic turbulent motions in accretion disks, especially for linearly stable fiows. A possible 
hydrodynamic origin of turbulent motions is important for several reasons. First, differential 
rotation is universally present in disks, whereas some disks or disk regions might not be ion- 
ized enough to support MHD phenomena, and a non-MHD source of turbulence must be at 
work there. For instance, protoplanetary disks are probably too resistive to support MHD 
turbulence (Fleming et al. 2000; Sano et al. 2000), but their observationally constrained 
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accretion rates imply the existence of turbulent transport. Secondly, the existence of self- 
consistent magnetized accretion/ejection structures seems to require a quasi equipartition 
of thermal and magnetic energy which, combined with the vertical stratification of these 
structures, may prevent the development of the magneto-rotational instability in a number 
of instances (Ferreira 1997; Casse and Ferreira 2000). 

The objectives of this paper are achieved through two different means. First, relevant 
pieces of information on the behavior of various types of shear flows which are available 
in the specialized fluid dynamical literature arc presented; from this material, it is argued 
that shearing sheet flows should be turbulent. Second, a phcnomcnological description and 
understanding of relevant turbulent properties of shear flows as they appear in the available 
experiments and numerical simulations is developed. In its most basic form, the phenomenol- 
ogy of hydrodynamic turbulent flows often relies on the concepts of Kolmogorov cascade and 
turbulent viscosity, and this approach is adopted here. Although the turbulent viscosity con- 
cept is of limited validity in complex situations (Tennckcs and Lumley 1972) and its use in 
the assessment of stability properties of turbulent flows has been rightly criticized (see, e.g., 
Terquem 2001; Hawley et al. 2001), it is well-known to provide accurate scalings of mean 
flow properties in simple shear flows such as channel or planar Couette flows (Tennekes and 
Lumley 1972; Lesieur 1987); consequently it has been widely used to parameterize turbulent 
transport in accretion disks. In this paper, some new and interesting consequences of this 
ansatz and their implications for turbulence in shear flows are pointed out. 

This paper is organized as follows. In the next section, some relevant flow conflgurations 
are introduced, along with the related forms of the Navier-Stokes equations; 1 also summarize 
relevant features of turbulence in these flows, as found in the literature, as this material 
has some direct bearing to the question of hydrodynamic turbulence in accretion disks, 
and is largely unknown to the astrophysical community; the reader who is not interested 
in factual details but only on their significance can jump directly to section 2.4, where 
this material is used to infer that "perfect" shearing sheet simulations should be turbulent. 
The most interesting findings of the present work are collected in section 3; after briefiy 
recalling the origin and rationale of the turbulent viscosity prescription, some of its previously 
unnoticed but important consequences are derived and used to interpret the behavior of 
the flows previously described, with special attention paid to Couette- Taylor flows, and 
to flow description in the shearing sheet approximation; in particular, a phenomenological 
explanation of the scaling of the Reynolds number with gap width in subcritical Couette- 
Taylor flows is devised. On the basis of this phenomenological understanding, the various 
reasons which might likely prevent the onset of turbulence in the simulations of Balbus et 
al. (1996) and Hawley et al. (1999) are discussed and the most critical one identifled. The 
final section summarizes the most relevant conclusions and discusses their consequences for 
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accretion disk theory and simulation, in particular on the magnitude of the Shakura-Sunyaev 
a parameter. 

The reader interested only in the new results of this paper and not in the background 
fluid mechanical information can focus on sections 2.4, 3.2 to 3.4, and 4. 

2. Turbulence in hydrodynamic shecir flows 

Hydrodynamic accretion disk mean flows are widely believed to be subcritical, i.e., the 
viscously relaxed laminar flow is linearly stable at all Reynolds numbers, at least locally. 
Furthermore, all experiments and numerical simulations of interest here pertain to linearly 
stable flows, and we are mostly interested in the local generation of turbulence. Therefore, 
I will focus on subcritical flows in this paper. 

The transition to turbulence is usually rather different in subcritical and supercritical 
flows. Supercritical flows undergo a cascade of precisely deflned bifurcations in parameter 
space, eventually leading to fully developed turbulence; these transitions are well docu- 
mented and reproduced numerically, e.g. for Couettc- Taylor flows (Andereck et al. 1986 and 
references therein; Marcus 1984a,b). Turbulence in subcritical flows, on the contrary, may 
abruptly be triggered, most probably by finite amplitude instabilities (Dauchot and Daviaud 
1994); also, the flow apparently evolves from highly intermittent to fully turbulent over a 
range of Reynolds numbers. 

Furthermore, shear flows can be either (wall-)bounded or free. The distinction refers 
to the limitation of the flow in the direction where the shear is applied (the transverse 
or shcarwise direction). This difference in boundary conditions inffuenccs some of their 
turbulent properties; indeed, free ffows are characterized by a single length-scale, the extent 
of the shear layer, whereas the distance to the wall introduces a second length scale in 
wall-bounded ffows. The inffuence of the other (streamwise and spanwise) boundaries is 
minimized inasmuch as their spacing exceeds the coherence length of the largest turbulent 
eddies, and as globally induced perturbations (such as Ekman circulation) are minimized by 
appropriate designs of the experimental setups. 

Shear ffows have been actively studied in the past decades, and their turbulent prop- 
erties are now characterized for a large variety of settings. In this section, I will brieffy 
present the subcritical ffows which have direct bearing to the question of hydrodynamic 
turbulence in accretion disks, namely, plane Couette and free shear ffows, either rotating 
or not, Couette- Taylor ffows, and Raylcigh-stable tidally driven shear ffows in the shearing 
sheet approximation. The ffrst two have been studied through both experiments and nu- 
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merical simulations. On the contrary, information on Rayleigh stable Taylor-Couette flows 
comes only from experiments. Finally, the shearing sheet approximation has been widely 
used as a local analytic model of local accretion disks, and has been implemented in the 
numerical work of Balbus, Hawley and coworkers quoted in the introduction (Balbus et al. 
1996; Hawley et al. 1999). For each of these flows, I characterize the geometry, the criti- 
cal parameters which are important for the question of the onset of turbulence, and I also 
give the governing dynamical equation (Navier-Stokes) in the form which is most suitable to 
establish comparisons between the various types of flows. 

The object of this section is to try to give an answer to the following question: if 
numerical simulations were perfect (i.e., not limited by questions of resolution, numerical 
instabilities etc) , would shear flows be turbulent in presence of the Coriohs force ? This is 
done in section 2.4, with the help of the material collected here. 



2.1. Plane Couette and free sheEir flows 

In spite of their conceptual simplicity, plane Couette flows are difficult to produce in 
actual experiments, which explains why some of their basic turbulent properties have only 
recently been characterized. The experimental setup is schematically represented on Fig. 1, 
along with a sketch of the turbulent mean flow proflle (see Tillmark and Alfredsson 1992 
for details). In practice the two walls are often made up of counter-moving (looped) inflnite 
belts. Similarly, free shear layers are produced by injecting fluid with different velocities on 
each side of a separating plate. The fluids come in contact at the end of the plate, and a 
turbulent layer develops and widens downstream (see Fig. 2). 

These flows are described by the Navier-Stokes equation in its simplest form, which 
reads 



+ v.Vv = + i/Av, (1) 

at p 

with obvious notations. The viscous terms are displayed in the incompressible form, as we 
are mostly concerned with subsonic turbulence. 

It is customary to deflne the Reynolds number of plane Couette flows based on the half- 
velocity difference (i.e. C/) and half-width distance (i.e. K) between the two walls. However, 
for the purpose of comparison with other setups, I shall deflne the Reynolds number as 
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Re = Wh/u, (2) 

i.e. based on the total velocity difference and distance between the two boundaries; the 
reader should bear in mind the resulting factor of 4 when comparing the figures quoted in 
this paper for Couette fiows with the literature. Prom the experiments of Tillmark and 

Alfredsson (1992), the minimal Reynolds for which turbulence is sustained is Re ~ 1500. 
The onset of turbulence in planar Couette flow has been successfully reproduced in numerical 
simulations (e.g., Bech et al. 1995 and references therein); a nonlinear mechanism for tapping 
the mean shear to sustain turbulence has even been identified (e.g., Jimenez and Moin 1991; 
Hamilton et al. 1995; Waleffe 1997). 

Rotating Couette and rotating free shear flows are produced by placing the experimental 
setup on a rotating platform. Such flows are very relevant to astrophysics, as they share a 
number of features with accretion disk flows in the shearing sheet approximation. For these 
fiows, the Navier-Stokes equation reads 



dv VP 

+ v.Vv = 2Q X V + Fin + i^Av, (3) 

at p 

where Fin. stands for the inertial force due to rotation. 

These rotating fiows are usually simulated by including only the Coriolis force term^ in 
the Navier-Stokes equation (Bech and Andersson 1996a; Komminaho et al. 1996; Bech and 
Andersson 1997), so that Eq. (3) reduces to 



(9v VP 

— + v.Vv = 2Q X v + i/Av. (4) 

at p 

Rotating Couette and rotating free shear fiows are characterized by the ratio S of the angular 
velocity of rotation to the shear 



2Q 

(5) 



d{v^)/dy' 



^The centrifugal term is not included on the basis that it results only in a redistribution of the equilibrium 
pressure. 
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where {vx) is the mean velocity profile; this number is akin to an inverse Rossby number, 
and measures the relative strength of the Coriolis and advection terms in the Navier-Stokes 
equation. A linear shear is destabilized by rotation when 

-1< 5 < 0, (6) 

and stabilized otherwise (see Tritton 1992 and references therein; see also section 2.4). The 
relevant regime for astrophysics is 5* < —1 (i.e., negative S and linearly stable linear shear^). 
No systematic exploration of the (i?e, S) parameter space has been performed; furthermore, I 
am not aware of any experimental investigation of rotating Couette flows for such (relatively) 
high values of S. However, this regime is explored in the set of free shear layers experiments 
of Bidokhti and Tritton (1992), who show that the flow remains turbulent^ (although linearly 
stable) down to 5* ~ —2 for Reynolds numbers^ ~ 4000 (see figures 14 and 16 of their paper). 
On the other hand, in the numerical simulations of anticyclonic {S < 0) rotating Couette 
flows of Bech and Andersson (1997) {Re ~ 5000) and Komminaho et al. (1996) {Re ~ 3000), 
turbulence is lost^ for ~ — 1 in the central part of the flow. This situation is similar 
to the one relating the simulations of Balbus et al. (1996) and Hawley et al. (1999) to the 
experimental data of Taylor (1936) and Wendt (1933) quoted in Richard and Zahn (1999); 
this analogy will be further discussed in section 2.4 and 3.4. 



^In relating rotating flows to shearing sheet ones, notice that the y axis identifies to the radial one, 
whereas the x and azimuthal directions are antiparallel. 

^By virtue of the Taylor-Proudman theorem, the flow should eventually become bidimensional but this 
happens only at higher values of l^l. 

^Notice that, as the turbulent shear layer widens downstream, Bidokhti and Tritton (1992) base their 
definition of the Reynolds number on the downstream distance x, which needs to be related to the layer 
width from which all Reynolds numbers quoted here are defined, and which is referred to as 25m in their 
paper. The two quantities can be related with the help of the various relations given in section 3 of their 
paper. This amounts to reducing the Reynolds numbers they quote by a factor ^ 7. Finally, the number 
given above corresponds to the most downstream point of measurement, where the flow should be closest to 
a developed (rather than developing) turbulent flow (incidentally, this is much farther downstream than the 
region where the pictures shown in the paper are taken). 

^Rotation in these numerical experiments is characterized by a global rotation number Rg = 2i}h/U 
rather than by the local rotation parameter S. In the central part of the profile, one usually has Rg > 0.2|S'| 
for fully turbulent flows, but it is difficult to precisely relate the relative level of rotation in these experiments 
to the critical limit between linearly stable and unstable rotating flows. 
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2.2. Couette- Taylor flows 

Couette-Taylor flows are produced from two concentric rotating cylinders. Most inves- 
tigations of this type have focused on the linearly unstable regime. The linearly stable one, 
which is more directly relevant to astrophysics, has only been explored by Taylor (1936), 
who maintained the inner cylinder at rest, and by Wendt (1933), who also reported results 
when the flow is close to marginal stability (i.e. linearly stable, but close to constant specific 
angular momentum) . The Reynolds number of these flows is defined as 



rAQAr 

Re = (7) 

where r is the mean of the two cylinder radii, and AQ and Ar are respectively the difference 
in angular velocity and the gap width of the two cylinders. Both investigations mentioned 
above did characterize the behavior of the minimal Reynolds number for well-developed 
turbulence to be maintained as a function of the cylinders relative gap width; this behavior 
is sketched on Fig. 3. Recently, a French team has undertaken an experimental investigation 
of flow profiles which are approximately keplerian in the mean, and found that turbulence 
was also maintained for Reynolds numbers of the order of a few thousand for a relative gap 
width of the order of a third (Richard, PhD thesis 2001). 

The minimal Reynolds number appearing in Fig. 3 is obtained by starting from an initially 
laminar flow, and progressively increasing the difference in angular velocity of the two cylin- 
ders (or only the outer cylinder angular velocity if the inner one is at rest). When starting 
from an initially turbulent flow and reversing the process, the loss of turbulence occurs for 
Reynolds numbers which can be significantly lower, but the fiow is then highly intermittent; 
it is reasonable to assume that the minimal Reynolds numbers of Fig. 3 are characteristic, 
albeit overestimated, values for well-developed turbulence (Richard, PhD thesis 2001). 

The two remarkable features of this minimal Reynolds number are a behavior which is 
similar to plane Couette fiows for Ar/r < 1/20, with Re ~ 2000, and a quadratic scaling 
[Re ~ Re*{Ar/rY with Re* ~ 6 x 10^] which is characteristic of rotation, as argued by 
Richard and Zahn (1999); these authors also show that in the same regime, the turbulent 
viscosity vt — I3r^\d^/dr\ with /3 ~ 10~^. A heuristic explanation of these features is 
presented in section 3. 

The Navier-Stokes equation for these flows is most meaningfully compared to that of 
other flows when substracting out the mean flow rotation JIq (i-©-, the average angular 
velocity of the two cylinders), as only differential rotation plays a role in the generation of 
turbulence. Deflning w = v — Jlo'^e,^, and (f) — 9 — Q^ot (so that w and are the velocity 
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and azimuthal coordinate in the rotating frame, respectively), the Navier-Stokes equation 
for w = {wr, Wrt>, w^) bccomes 



— + w.V'w + 2Q X w H —es = h Br + i/Aw, (8) 

at r r \ p J 

where w.V'w = (w.Vu'r)er + (w.Vw0)e<p + ((w.Vwz)^^. For future reference, I refer to 
the terms w'^/r and 2wrWo/r as "geometric terms", as they arise because of the cyhndrical 
geometry^. 

For these flows, the rotation parameter defined in Eq. (5) reads 



rdfl/dr g' 

where q = —{r/Q){dQ/dr) is the parameter defined by Balbus et al. (1996) to characterize 
rotation profiles. The flow is stable according to Rayleigh's criterion when q < 2, i.e. when 
S < —1, quite similarly to rotating Couette and free flows, although the processes through 
which instabihty occurs are different. Note also that Eqs. (4) and (8) differ only through the 
geometric and centrifugal terms. 

The fact that the minimum Reynolds number for developed turbulence is identical in 
plane Couette and Couette- Taylor flows with Ar/r < 1/20 and the inner cyhnder at rest 
can be understood in the following way. First, the advection term (which is the source of the 
turbulence cascade as indicated by the very existence of the Reynolds number) dominates 
over the geometric terms when rAQ/Ar 3> rAQ/r, i.e. Ar/r <S 1. Second, Afl = fl (one 
cylinder being at rest), so that the Coriolis term is also very small compared to the advection 
term, and Eq. (8) nearly reduces to Eq. (1). Note furthermore that the Coriolis force does 
not appear to significantly affect the minimal Reynolds number for the onset of turbulence 
for the values of q of interest here (i.e., g^ltog~l-2), both in the limiting plane 
Couette regime and in the rotation regime, as exemplified by the data of Wendt (1933) for 
nearly neutral flows, which follow the same law for the minimal Reynolds number, down to 
the plane Couette limit. 



^Such terms also arise in principle from the viscous term, but they are inessential to the argument 
developed in this paper. 
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2.3. Shearing sheet 

Accretion disk flows in the shearing sheet approximation are closely related to the 
Couette- Taylor flows previously described. They differ in only three respects. 

First, the mean angular velocity profile (Q) in Couette-Taylor flows is imposed by the 
boundary conditions, and by the condition of stationarity of the mean viscous or turbulent 
transport of angular momentum, from the azimuthal momentum equation (with the walls act- 
ing as source and sink of angular momentum). The radial momentum equation then imposes 
the mean radial pressure profile (P), and the resulting tidal force term (— V(P)/(p) +r(f2)^). 
On the contrary, in accretion disks, the mean radial angular velocity profile mostly results 
from the gravitational attraction of the central body, which imposes a nearly keplerian pro- 
file in cold disks, but the disk is never globally stationary, due to viscous/turbulent trans- 
port (nevertheless, an approximate stationarity is nearly achieved locally on the dynamical 
timescales of interest for the onset of turbulence). Therefore, in keplerian disks in the shear- 
ing sheet approximation, the tidal force term {—g + r(0)^) is the source of the (keplerian) 
angular velocity profile and not its consequence. Furthermore, one usually neglects the radial 
pressure gradient locally, and assumes that the gravitational force has cylindrical (and not 
spherical) symmetry for simplicity, as cold disks are thin. 

Secondly, a local approximation is performed, by restricting consideration to a radial 
box of width Ar <^ r; one also usually assumes that the height of the box is comparable 
to its width. Under these assumptions, one neglects the geometric terms in Eq. (8), and 
describes the fiow in local cartesian coordinates {x <->■ r, y <->■ r(f) where is the azimuthal 
coordinate in the rotating frame introduced for Couette-Taylor fiows). One also linearizes 
the angular velocity profile. 

Finally, this local approximation and the resulting change of geometry from cylindrical 
to cartesian (except for the Coriolis force term which is kept) allows one, in numerical 
simulations, to adopt a particular form of the periodic boundary condition in the radial 
direction, in which the fiuid quantities on the radial boundaries are longitudinally displaced 
all the time with the mean angular velocity difference during a time step before the periodic 
boundary condition is applied (see Hawley et al. 1995 for details on this procedure). 

With these prescriptions (aside from the boundary conditions) , the Navier-Stokes equa- 
tion, in the shearing sheet approximation and in the rotating frame, becomes 



— w. Vw 2Q X w = h 2qn'^xer + i^Aw, 

at p 



(10) 
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where x — r — tq and tq is the position of the center of the shearing sheet box. The term 
2qQ'^x represents the tidal force (difference of the gravitational and inertial force); q is the 
parameter introduced in the previous subsection and measures the steepness of the rotation 
profile. Note that the pressure term contains only fiuctuations related to the presence of 
turbulence, which is not the case in Couette- Taylor fiows. It is interesting to note that this 
equation shares features with both Eqs. (4) and (8); in particular, linear stability is ensured 
for < 2, i.e. < — 1 for the laminar linear profile. This makes the loss of turbulence in 
the simulations of Balbus et al. (1996) and Hawley et al. (1999), for values of q smaller than 
2 by a few percents only, all the more intriguing. 



2.4. Shearing sheet, rotating Couette flows, and turbulence 

In fact, all available pieces of evidence strongly suggest that numerical simulations of 
rotating Couette fiows and of tidally driven sheared motions in the shearing sheet limit 
should display turbulence, as I argue now. 

First, plane Couette flows, rotating Couette flows, and tidally driven shearing sheet flows 
have similar linear stability properties. For all three types of flows, the viscously relaxed 
laminar solution is a simple linear shear, which is always linearly stable for the plane Couette 
fiow^, and stable for the other two fiows once 5" < — 1 (which is the only case of interest 
here). Plane Couette fiows are subject to finite amphtude instabihties (see, e.g., Lerner 
and Knobloch 1998, DubruUe and Zahn 1991, and references therein). The same is true of 
rotating Couette flows (Johnson 1963), and of shearing sheet fiows (Dubrulle 1993). As finite 
amplitude instabilities are considered to trigger the turbulence seen both in experimental 
and numerical investigations of plane Couette flow, one would expect the same to be true of 
the other two flows. 

Secondly, let us reexamine the differences between rotating Couette flows and the shear- 
ing sheet flows with the other flows discussed previously. They amount to differences in 
boundary conditions, of mean force terms, and of geometry. 

The shearing sheet boundary conditions are in a way intermediate between rigid and 
free boundary conditions, as they imply that the mean flow obeys rigid boundary condi- 
tions, whereas the fluctuating part obeys periodic boundary conditions; rotating Couette 
flow simulations are usually performed with rigid boundary conditions. On the one hand 



consider unbounded flows in this discussion, as instabilities due to the boundary in viscous fluids are 
not relevant in astrophysics. 
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Couette- Taylor flows implement rigid boundary conditions. Although in real experiments, 
the vibrations of the boundary play some role in triggering turbulent motions, there is lit- 
tle doubt that in these experiments, turbulence is self-sustained. On the other hand, from 
the experiments of Bidokhti and Tritton (1992) the rotating shear flows with free boundary 
conditions are turbulent, although by construction no mean steady state can be reached in 
these systems. Therefore, it seems unlikely that boundary conditions play an important role 
in the presence or absence of turbulence in numerical experiments. 

In the shearing sheet approximation, the mean shear is imposed by the tidal force term; 
in rotating Couette simulations, it results from the boundary conditions. In Couette- Taylor 
flows, the boundary conditions do not only produce the shear, but also generate a mean 
radial pressure gradient. Note however that the term —d{P)/dr/{p) + r(f2)^ of Eq. (8) is 
similar in function to the term 2qQ,^x in Eq. (10). Furthermore, the mean pressure gradient in 
Couette- Taylor experiments is radial, whereas it is longitudinal (streamwise) in the rotating 
free shear layer experiments of Bidokhti and Tritton (1992). This suggests that neither 
large scale mean pressure gradients, nor tidal terms, make any significant difference on the 
question of the onset of turbulence in the various flows considered here, especially that all 
gradient terms get out of the way in incompressible flows (they disappear from the vorticity 
equation) . 

Finally, 1 will show in the next section that the main effect of the geometry (which 
enters through the geometric terms in Couette- Taylor flows) is to change the conditions of 
onset of turbulence, but this does not affect the occurrence of turbulence in itself. 

Although such arguments do not exclude more complex possibilities (as, e.g., that tur- 
bulence might be impeded in shearing sheet flows by a combination of these factors instead 
of only one of them), this strongly indicates that rotating Couette flows and shearing sheet 
ones should be turbulent, suggesting that the absence of turbulence in all the published 
simulations of this kind stems from limitations in the numerics involved. This last point is 
addressed in the next section. 



3. Phenomenology of subcritical turbulence 

The purpose of this section is to point out important features of turbulence in sheared 
flows, through a phenomenological model developed in section 3.2. The consequences of 
this model are used in section 3.4 to identify the potential limitations in the numerics just 
mentioned. 
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3.1. Turbulent viscosity and the Kolmogorov prescription 

In a picture where the fluctuating turbulent scales can be separated from the more 
regular large ones, it is meaningful to write down an equation for both the mean (X) and 
fluctuating 6X parts of any quantity X. In particular, the evolution of the mean velocity 
reads 



^ + V • (SySv) = + uA{y), (11) 

ot p 

where possible geometric and/or inertial terms have been omitted for simplicity, as well as 
the effect of compressibility. In the simple shear configurations of interest here, only the 
(SvySvx) (or (SvrSv^)) part of the Reynolds stress tensor is relevant for radial turbulent 
transport. 

By describing turbulent fiuctuations with a characteristic coherence scale Im and velocity 
amplitude vm, Prandtl (1925) argued that 



{SVySV^ 



dy 



(12) 



with 



~ ImVm- (13) 

Note that in cylindrical geometry {5vr5v^) ~ Ufrd^D,) /dr. 

The reasoning behind this formulation is similar to the one relating the usual molecular 
viscosity to the molecular mean free path and velocity dispersion (i.e., turbulent transport 
occurs over a "mean free path" Im with "velocity dispersion" vm)', Eq. (12) can also be derived 
from more rigorous multi-scale expansion techniques. In a Kolmogorov cascade picture, Im 
is the energy injection scale (and characterizes the coherence length of the largest eddies of 
the cascade), and vm the amplitude of the velocity fluctuations at this scale, as the velocity 
amphtude decreases with decreasing scale in a Kolmogorov spectrum. However, this does 
not mean that larger fluctuating scales are not present in the flow, nor that they have no 
influence in the development of turbulence; it just implies that they dominate neither the 
energy spectrum nor the turbulent transport. 



An important feature of the turbulent viscosity prescription is that the rate of energy 
dissipation e — which is also the rate of energy transfer in eddy-scale (Fourier) space — is 
simply given by 



^ ~ Vm/Im 




{d{vx) I dyf' (cartesian) , 
Vt {rd{'D,) I dr^ (cylindrical) , 



(14) 



as can be shown most directly by deriving the relevant macroscopic energy equation. Eqs. (12) 
and (14) imply in particular that the characteristic frequency of turbulent motions is the 
shear frequency, i.e. 



This reflects the fact that an externally imposed shear locally possesses no characteristic 
scale (besides the scale of the flow), but only a characteristic frequency, so that shear turbu- 
lence can only couple efficiently to the shear if its characteristic frequency or coherence time 
(at the energy injection scale imposed by the mechanism which drives turbulence) matches 
the shear frequency®. 

The turbulent viscosity description has been applied to a wide variety of setups to 
describe the mean properties of turbulent flows, both in the the vicinity of walls and in the 
main part of either free or bounded flows (e.g., Tennekes and Lumley 1972; Lesieur 1987). 

3.2. Turbulence scales : phenomenological model and orders of magnitude 

My primary purpose here is to point out some interesting consequences of the turbulent 
viscosity prescription. Indeed, one expects that a flow undergoes a transition to turbTilcnce 
when the turbulent transport becomes more efficient than the laminar one for subcritical 
flows. This implies that 



*As the coherence time of smahcr scale eddies is shorter, they are less or little affected by the shear. As a 
consequence, in a first approximation, the turbulence is more or less isotropic at scales < Im, and anisotropy 
is ignored in the whole argument. 




(cartesian) , 
(cylindrical) . 



(15) 
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i/j > 1/ when Re > Rem, (16) 

where Re^ stands for the minimum Reynolds numbers for the onset of turbulence discussed 
in the previous section. For example, note that for Couette- Taylor flows, from the data of 
Wendt (1933) and Taylor (1936) z/j/i/ ~ l3Re* ~ 6, where (3 and Re* are the quantities 
introduced in the previous section in the discussion of these flows; furthermore, when the 
minimal Reynolds number is searched for by decreasing the velocity difference between the 
cylinder from an initially turbulent state instead of increasing it from an initially laminar 
one, the ratio i^t/i^ is sensibly closer to unity (Richard, PhD thesis 2001). 

Away from boundary layers (if any) , the only scales which are relevant for characterizing 
the shear gradient are the typical size of the shear flow. Ay (resp. Ar) in cartesian (resp. 
cylindrical) geometry, and the typical shear amplitude over this scale, Af^ (resp. rAQ). 
Combining Eqs. (12), (14) and (16) then yields, for the bulk of the turbulent flow 



I'M 



(Ay/Rel^^ (cartesian), ^^^^ 
l^Ar/Rell"^ (cylindrical). 



and 



jAvcc/ReJ^ (cartesian), 
[rAQ/Rell^ (cylindrical). 

I wish to stress that Eqs. (17) and (18) do not imply that turbulence is a global rather 
than local phenomenon. On the contrary, Eqs. (14) and (15) relate Im and vm to local char- 
acteristics of the mean flow. Note that these relations justify (at least for subcritical flows) 
the separation of scales between the mean large scale flow and the fluctuating small scale 
one which is assumed in the turbulent viscosity description, because Rem usually exceeds a 
few thousands. 

These relations have a direct physical interpretation. Consider for example two planar 
Couette flows with identical shear rates, and with wall spacing Ay and relative velocity Av^ 
which differ by a given ratio. Obviously, the scaling with Ay and Avx is a natural conse- 
quence of the scaling similarity between flows which are otherwise identical. On the other 
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hand, consider different flows, with identical shear rates, but different minimal Reynolds 
numbers (e.g., plane Couette and Couette- Taylor flows with appropriate parameters). A 
larger minimal Reynolds number is a sign of a greater difficulty to trigger turbulence, i.e. an 
increased difficulty for turbulent transport to dominate over the viscous one, and therefore 
is a sign of a smaller scale turbulence, due to the physical picture underlying the turbulent 
viscosity prescription (i.e., the transport occurs over a smaller "mean free path" Im, and 
correlatively with a smaller "random velocity" vm due to the assumption of identical shear 
rate between the two different ffows). 

From these relations, one can easily check that, at the minimum Reynolds number, 
the advection term, which dominates scale coupling and is the primary cause of the incrtial 
turbulent spectrum, is comparable to the dissipation term, at the turbulent transport scale. 
As a consequence, the turbulence possesses little or no inertial domain at its threshold. 
Furthermore, as long as there is no change in the turbulence generating process, increasing 
the Reynolds number can only result in lowering the dissipation scale with respect to Im, 
and therefore in the progressive build up of an inertial spectrum (e.g., imagine one does this 
by reducing the viscosity while maintaining the large scale structure of the ffow unchanged) . 

It is important to notice that the estimates of Eqs. (17) and (18) remain vahd for 
Reynolds numbers larger than the turbulence threshold, as long as the turbulence generating 
process is unchanged. The predictions of the scaling proposed here are well supported by 
the available empirical and numerical evidence, as shown in Appendix A. 

3.3. Consequences: Couette- Taylor flows 

Eqs. (17) and (18) have particularly interesting consequences for the understanding of 
turbulence in Coucttc-Taylor flows. For dcfinitcness, I will flrst focus on flows where the 
inner cylinder is at rest. As argued at the end of section 2.2, for r ^ Ar, the Navicr-Stokes 
equation for Couctte-Taylor flows [Eq. (8)] then reduces to the Navier-Stokcs equation for 
planar Couette flows [Eq. (1)] and the minimal Reynolds number is constant. However, when 
Ar — >■ r, the geometric terms 0{w'^/r) ~ {rAfl)'^/r become comparable to the advection one 
on scale Ar. Furthermore, if, at some radial location r in the flow, Rcm remained constant 
when Ar ^ r, Eq. (17) would imply that Im could become arbitrarily larger than r, which 
makes little sense. In fact, one expects that Im oc r once Ar/r exceeds some critical ratio 
Ac (which for the time being is expected to be of order unity), for two reasons: flrst, the 
geometric terms introduce a limiting scale (the radius r), which must be accounted for by 
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the turbulent viscosity description^; second this prescription for Im is necessary to satisfy 
the requirement that luVw > w'^/r at the largest scale of the inertial spectrum (in order to 
maintain such a spectrum). Consequently, let us assume that 



Im ~ jr, (19) 

when Ar/r > Ac, and where 7 is a constant to be determined later. The argument presented 
here suggests that Ac ~ 1 whereas the data imply that it is significantly smaller than unity 
(see below); as for the large values of the minimal Reynolds numbers for turbulence (which 
one would also naively expect to be of order unity), this originates in the (still unknown) 
mechanism which sustains turbulence. 

Eqs. (17) and (19) must be satisfied simultaneously, and this is possible only if 
depends on the relative gap width: 



Rem-\(—)\ (20) 



which explains the behavior seen on Fig. (3). Equivalently^*^, r'^ {(El / dr) / v > 1/7^; this shows 
that, as soon as Ar > Ac, the width of the fiow does not influence the onset of turbulence, 
which becomes a purely local phenomenon. 

The velocity fluctuation amplitude now reads 



'"M ~ 1-^ - ir — , (21) 

i.e., it is proportional to the local shear rate. As a consequence, the turbulent viscosity 
becomes 



^This reasoning is somewhat similar to the one which imposes that oc y in the vicinity of the wall in 
Couette or channel flows, and which has lead to the derivation of the well-known "law of the wall" , describing 
the mean structure of turbulent flows close to the wall (e.g., Landau and Lifshitz 1987; Tennekes and Lumley 
1972; Lesieur 1987). 

^•^Remember that Ar and Afi have been introduced in Eqs. (17) and (18) to represent local gradients in 
order of magnitude. 
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(22) 

A similar relation has also been proposed by Richard and Zahn (1999) directly from exper- 
imental torque data. Note that the reasoning leading to Eq. (22) implicitly assumes that 
the flow compression plays little role, so that this result may not necessarily apply to super- 
sonic turbulence. In order for Eqs. (20) and (22) to faithfully account for the properties'^ of 
Couette- Taylor flows described in section 2.2, one needs to tie up a few loose ends: 



• Because the Coriohs force does not seem to affect the minimal Reynolds number of 
turbulence (see the closing comment of section 2.2), the argument above must apply 
to any value of the q < 2 (the parameter introduced in section 2.3 to characterize 
the local rotation proflle), and not only to situations with the inner cylinder at rest; 
however, for g ~ 1, the geometric term is always comparable to the advection term, 
and the argument is less transparent. 

• The gap relative width in Fig. 3 is measured with respect to the mean radius of the 
rotating cylinders, whereas a local value is used above. However, the relative gap widths 
shown in this flgure are all sufficiently smaller than unity to make the difference between 
the two quantities negligible in the scaling argument developed here. Incidentally, 
this shows again that turbulent properties are local; e.g., turbulent eddies become 
larger when one moves outwards in a sufficiently wide cylindrical system. Correlatively, 
Eq. (19) follows also directly from the fact that r is the only available local scale. 

• The relations derived above imply that 7^ = /? = Re*~^, but the values quoted in 
section 2.2 for the last two quantities differ by a factor of 6. However, it was also 
pointed out there that the value of Re* ~ 6 x 10^ is overestimated, because it leads 
to a ratio Ut/u which is too large due to the particular experimental protocol adopted 
by Taylor (1936) and Wendt (1933). Also, recent (still unpublished) experiments on 
Couette- Taylor flow in the "keplerian" regime {q = 3/2) exhibit sustained turbulence 
for Reynolds numbers smaller than the limit of Fig. 2, but with a different experimental 
protocol (Richard, PhD thesis 2001). Therefore, it is reasonable to assume that /3 is a 
much better measure of 7^ than Re*; this assumption is made in the remainder of this 
paper, and (3 is used everywhere instead of 7^. 



^^Most notably, the scaling of the Reynolds number with (Ar/r)^ in the rotation regime, and the near 
coincidence between /3 and Re*~^. 
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The critical value of the relative gap width Ac which separates the planar regime from 
the rotating one obtains when the values of Im in both regimes are equal. This yields 



where Rcp ~ 2000 is the minimal Reynolds number in the planar limit. This gives Ac ~ 1/7, 
which is somewhat larger than the value of 1/20 shown on Fig. 3 (but closer to the uneducated 
guess Ac ~ 1), because of the reduction adopted above of the value of Re*. 

Note also that Im ~ r/300. One might wonder why such a small length scale arises, 
whereas one would naively expect ~ r on dimensional grounds. However, the same 
dimensional type of argument would also predict that turbulence sets in for Re > 1, which 
is strongly violated by the empirical evidence. The two facts have the same physical origin: 
the (as yet not understood) mechanism which sustains turbulence. 

Two other explanations of the behavior of the Reynolds number with relative gap width 
have previously been proposed in the literature. Zeldovich (1981) assumed that turbulence 
in these Couette-Taylor flows is controlled by a competition between the cpicyclic (stabi- 
lizing) frequency and the shear rate which is the source of the turbulent motions; however 
his findings are inconsistent with some of the data (sec the discussion of this point in the 
appendix of Richard and Zahn 1999). Dubrulle (1993) looks for an explanation in terms of 
finite amplitude instabilities in the WKB approximation, but this is incompatible with the 
fact that the scale r plays a key role in the problem. 

I conclude this section by pointing out that the Coriolis force appears nowhere in the 
arguments presented in this section, which suggests that it plays little role in the development 
of turbulence in subcritical Couette-Taylor fiows, at least for g ~ 1 — 2. Indeed, in opposition 
to the inertial (geometric) terms, the Coriolis force does not single out any length scale. In 
particular, the ratio of the advection term (~ rnVw) to the Coriohs one (~ wfl) in Eq. (8) 
is ~ 1 both at scale Im and at scale Ar for the values |g| ~ 1 — 2 of interest here, and 
increases with decreasing scale in a Kolmogorov cascade picture. However, it does play a 
role in the loss of turbulence in simulated rotating flows, but this apparent paradox cannot 
be investigated in the framework of the order of magnitude arguments developed in this 
section. The next section is devoted to a discussion of this point. 




(23) 
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3.4. The role of the CorioHs force: beyond orders of magnitude 

The question I want to address here is the following : why is turbulence lost in numerical 
simulations of sheared flows when even a small amount of rotation is added (and the result- 
ing flow remaining linearly stable) — in particular for Couette flows and disk flows in the 
shearing sheet approximation — whereas in experiments as different as the Couette- Taylor 
and rotating free shear layers, it is maintained (in the same conditions of linear stability). 
The forms of the Navier-Stokes equation for these flows given in section 2 strongly suggests 
that this is a consequence of the Coriolis force, as this is the only new force term which is 
taken into account when rotation is added to free shear layers and planar Couette flows in 
these simulations. 

More specifically, clues to the role of the Coriolis force can be found by inspecting 
the behavior of plane Couette and free shear flows with and without rotation. The most 
noticeable and important feature is that, in numerical simulations of rotating Couette flows, 
even a small Coriolis term is able to suppress the very largest scales of the turbulent motions. 
This is particularly obvious, e.g., when comparing Figs. 7 and 20 of Komminaho et al. 

(1996) , which shows that the very large scales which develop in the streamwise direction in 
turbulent flows break up for rotation numbers as small as a few percents. This feature is 
quite understandable on the basis of the velocity spectra shown in Bech et al. (1995), which 
imply that kv^ is most likely sensibly smaller than Q at scales larger than Im- A similar 
feature can also be indirectly found in Bidokhti and Tritton (1992) (see their Fig. 11 and 
14, to be combined with their Fig. 16), who show that the Reynolds stress tensor magnitude 
decreases by a factor^^ of at least 10 when the parameter S introduced in Eq. (5) varies from 
to < — 1; this suggests that the size of the largest turbulent scales in these flows is also 
substantially reduced under the action of the Coriolis force^^. This indicates that, although 
Eq. (17) always provides reliable orders of magnitude for Im, it underestimates the relevant 
eddy scale by a factor ~ 3 for non-rotating flows (see Appendix A), while overestimating 
it by at least the same factor once rotation is introduced. As a consequence the loss of 
turbulence in the numerical simulations of rotating Couette flows of Bech and Andersson 

(1997) and Komminaho et al. (1996) is clearly an effect of the limited small scale resolution 
due to the large box sizes (especially in the streamwise direction) adopted in these works: 



^^Thc noise in the data at large rotation number docs not permit a very precise estimate of this reduction 
factor, but the value of 10 quoted here seems a bare minimum. 

^^Note that in Bidokhti and Tritton (1992), as pointed out by the authors themselves, turbulence is not 
lost, and the flow remains three-dimensional, although the velocity fluctuations anisotropy is clearly affected 
by rotation. 
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the smallest available scales do not allow these authors to account for the inertial part of 
the energy spectrum, while all the larger scales are wiped out by the Coriolis force. It is 
more than likely that, in these simulations, the Coriolis force kills the large scale mechanism 
which has been identified to sustain turbulence in plane Coucttc flow simulations (Jimenez 
and Moin 1991; Waleffe 1997; Hamilton et al. 1995; see section 2.1). The fact that free 
rotating layers and Couette- Taylor flows remain turbulent at larger levels of rotation than 
the ones to which turbulence is lost in these simulations implies that a different mechanism 
for sustaining turbulence is at work in these flows, and that it operates at scales comparable 
to, but apparently smaller than, the estimate of Eq. (17). This other mechanism has not yet 
been found in numerical simulations. It would be interesting to know whether this change 
of mechanism is related to the fact that the Coriolis force apparently selects the direction of 
instability of flnite amplitude defects (Johnson 1963). 

The same line of argument applies to the shearing sheet simulations of Balbus et al. 
(1996) and Hawley et al. (1999). Indeed, the effective Reynolds number of these simulations 
is not an issue, as the code used by Hawley et al. (1995) and Hawley et al. (1999) is able to 
flnd turbulence — or at least the large scale mechanism already alluded to — in non-rotating 
Couette flows, and this happens only for Reynolds numbers larger than at least 1500. Also, 
the argument developed in section 3.2 shows that the Coriolis force by itself should not change 
the minimal Reynolds number for the onset of turbulence, an inference confirmed by the fact 
that turbulence is seen developing in the rotating free shear layer of experiments of Bidokhti 
and Tritton (1992) for roughly comparable Reynolds numbers. Under the assumption (cf 
the arguments developed above) that Eq. (17) provides an estimate for the largest turbulent 
scale which is overestimated by a factor of at least 3 in the presence of a Coriolis force term, 
one obtains Im ^ Ay/lOO, with (possibly much) smaller values more than likely. This is 
most probably too close to the largest resolution achieved in the shearing sheet simulations 
(the number of zones in any direction being no larger than 250), especially when artificial 
viscosity is taken into account, for turbulence to show up in these simulations. 

To conclude this section, it is worth noting that some other numerical questions must 
be considered to find turbulence in these systems. First, it is well known from experiments 
with subcritical fiows that the way perturbations of the fiow are designed has an infiuence 
on the appearance of turbulence. This suggests that some care must be exercised in the 
choice of the initial conditions in numerical experiments; in particular, it might be useful 
to ensure that at least some condition of finite amplitude instability is satisfied in this 
choice. Secondly, the role of the choice of the Courant number is not completely obvious, 
even in situations where the CFL condition is not violated. For example, in a series of 
yet unpublished simulations of linearly stable Couette- Taylor fiows performed with the Zeus 
code in collaboration with David Clarke, we did initially find that turbulence would set in for 
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flows which are "not too far" from planar Couette flows (including some roughly keplerian 
flows), but it would eventually disappear when reducing the maximal allowed time-step, 
although the CFL condition was satisfied in all runs. The reason of this behavior is not 
yet completely elucidated, but it appears to have some direct connection to the question 
of resolution just discussed^^. In any case, the disappearance of turbulence in numerical 
simulations of Couette- Taylor flows which are experimentally known to be turbulent is a 
serious cause of worry on the reliability of the conclusions drawn from the published shearing 
sheet numerical experiments. 



4. Summciry and Eistrophysical implications 

Eqs. (17) and (18) (section 3.2) along with their consequences constitute the central 
findings of this paper. They express the natural length and velocity scales which are in- 
volved in turbulent transport in subcritical flows in terms of the local mean characteristics 
of the flow, and result from the constraint that the turbulent transport dominates over the 
viscous one in the framework of the turbulent viscosity description. The scaling and orders 
of magnitudes implied by these relations is supported by the available experimental and 
numerical evidence (see the appendix and the beginning of section 3.4). 

These scaling relations have two important consequences. First they provide an ex- 
planation for the minimum Reynolds number dependence on the relative gap width in 
Couette- Taylor experiments, displayed on Fig. 3 (section 3.3). A theoretical explanation 
of this behavior has long been sought for, but none has satisfyingly been proposed yet; the 
phenomenological one presented here has the advantage of connecting apparently unrelated 
features, to be consistent with all the experimental constraints, and to point out the direction 
in which such a theoretical explanation might be looked for. Incidentally, the existence of 
this phenomenological explanations strengthens the validity of these scaling laws. Secondly, 
the comparison of the various flows presented in section 2 implies that disk flows described 
in the framework of the shearing sheet approximation should be turbulent (section 2.4), and 
the scaling relations strongly suggest that the absence of turbulence in the available shearing 
sheet numerical simulations is due to a lack of resolution (section 3.4). This follows because 
the Coriolis force destroys large scale fluctuations, thereby affecting in a major way the non- 
linear mechanism through which turbulence is maintained. At present, this mechanism is 
not understood, except, to some extent, for plane Couette flows. 



^"^A counterintuitive dependence of hydrodynamic simulations on the Courant number is also visible on 
Fig. 1 of Porter and Woodward (1994). 
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Understanding to which extent these results are helpful in characterizing and quantifying 
turbulent transport in accretion disks is an important issue. Three factors at least must be 
accounted for: the magnitude of the disk pressure, the vertical scale height, and the presence 
of a magnetic field; these factors are not macroscopically independent, but relate differently 
to the onset of turbulence. 

The disk pressure affects the problem in two a priori different ways: first, the turbulent 
transport picture presented in this paper requires the underlying turbulence to be subsonic 
(see also Hure et al. 2001); second, turbulent velocity fluctuations require a force to produce 
them, and only the pressure force is available to this purpose in the hydrodynamical case, 
independently of the details of the underlying mechanism which sustains this turbulence. 
The first constraint is easily quantified: turbulent motions are subsonic if vm/cs < 1 (c^ is 
the sound speed); in accretion disks, Cg ~ and from Eq. (15), this implies that H >Im 
{H is the disk scale height). To quantify the second constraint, note that a given fluctuating 
blob of size Im undergoes a velocity change Su ~ luirdQ./ dr) ~ Im^ over a time-scale 
~ vm/Im ~ ^) because the coupling to the shear is the source of turbulent motions at the 
largest scales; the largest pressure variation at any scale is 5P/ p ~ c^, and requiring that the 
resulting pressure force at scale Im is able to account for the turbulent velocity fluctuations 
at this scale requires^^ again H > Im- 

The turbulent scales (~ Im) are connected to the mean flow scales through the mech- 
anism which sustains turbulence. In an accretion disk, only two such mean flow scales are 
available locally: H and r. The role of r has already been discussed; the role of the vertical 
scale height depends on the anisotropy of the mechanism which sustains turbulence. In the 
absence of constraint on the nature of this mechanism for rotating shear flows, I will examine 
in turn two limiting assumptions: 

• This process is "isotropic", i.e. the scales it requires to operate are roughly identical 
in all directions — shearwise, streamwise and spanwise (this is the case for example 
of the nonlinear mechanism mentioned in section 2.1 for non-rotating plane Couette 
flows). In this case, the elementary box in which this mechanism operates must be 
of size H, which implies in particular that Ar ~ iJ in all the relations used in the 
previous sections of this paper. If H/r < Ac [cf Eq. (23)], as expected in most disk 
models, the rotation regime of Couette- Taylor flows is irrelevant; instead, the shearing 
sheet approximation appUes. As argued at the end of section 2.2, the Coriolis force is 
not expected to affect sensibly the minimal Reynolds number of turbulence, so that 



This argument ignores the possibiUty of supersonic turbulence. 
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ut ~ {rAfl)H/Rep ~ 10 ^CgH, and the Shakura-Sunyaev parameter a ~ 10 ^. Note 
in this case that the constraint H >Im is always satisfied. 

• The process is not sensitive to the vertical scale height except through the pressure 
requirement described above. As a consequence, as long as H > ~ 3.10~^r (which 
is hkely to be satisfied in accretion disks), the Couette-Taylor rotation regime applies, 
and Ut ~ Pr^dQ/dr ~ /3Qr^, so that the Shakura-Sunyaev parameter a ~ P{r/Hy lies 
in the range 10-=^ - 10"^ If H < p^/'^r, the turbulence is supersonic. Note however 
that the extra energy dissipation taking place in shocks makes a supersonic turbulence 
more difficult to maintain, and the disk might heat up until ~ /jv/ is satisfied again, 
or Im might decrease, i.e. the turbulence maintaining process might be afi'cctcd and 
increased, or the limiting case considered here does not apply^^. This makes the 
relevance of supersonic turbulence to accretion disk theory unclear. 

The conclusion of this brief discussion is that in the Shakura-Sunyaev parameterization 
of the turbulent viscosity in hydrodynamic disks, either a ~ 10~^, or a ~ 10~^(r/iJ)^, 
depending on the unknown characteristics of the mechanism which sustains turbulence. In 
principle, one should also check that Re > Rem', as z/ ~ /c (c ~ Cg is the velocity dispersion 
and / the mean free path), this translates into H/l > 10^ in the first case above, and /HI > 
j3^^ ~ 10^ in the second, but both requirements are most probably satisfied everywhere in 
astrophysical accretion disks. 

It is unclear how the presence of a magnetic field can modify hydrodynamic shear turbu- 
lence. In particular, even a dynamically non dominant field can easily affect the mechanism 
of generation of turbulence, and therefore significantly modify the efficiency of the turbulent 
viscosity transport, on top of adding a turbulent resistivity, even if the MHD fiow remains 
linearly stable^^. Reversely, the possible occurrence of hydrodynamic shear turbulence can 
possibly affect in a major way our present understanding of MHD transport and dynamo 
processes in accretion disks, which mostly relies on the physics of the nonlinear development 
of the magneto-rotational instability. Clarifying these questions is of primary importance 
for accretion disk theory. 



"^^If the Reynolds number is large enough, the disk must be turbulent; this follows by considering a narrow 
enough disk portion so that H exceeds its width, and at least one of the regimes of the previous sections 
does apply, inasmuch as boundary conditions arc not essential to the onset of turbulence, as argued above 

^^This can happen, e.g., if the disk scale height is small enough as not to let any magneto-rotational mode 
become unstable, which is easily realized in disks with a near equipartition between thermal and magnetic 
energies. 
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To conclude this paper, let me point out that there is one example of keplerian disk 
which has been observed with a great luxury of details, and which is not turbulent, namely 
Saturn's rings. However, the requirements discussed above fail on several accounts in ring 
systems, because both the particle size d and mean three path / are comparable to H . For 
example, in the first limiting case discussed above, the ring is necessarily laminar, while in 
the second, because H > Im, the granularity of the system makes scales < Im inaccessible 
to the fluid description; the same argument makes supersonic turbulence most probably 
irrelevant to ring systems. 



A. Evidence for the proposed turbulent scaling 

Because Rcm is at least of the order of 10'^, the order of magnitude estimates of Eqs. (17) 
and (18) are sensibly smaller than the mean flow length and velocity scales to which Im and 
vm are usually assumed to be comparable. Nevertheless, they are in good order of magnitude 
agreement with the available evidence. Consider, for example, the simulation of Couette 
flow reported in Bech et al. (1995), and further exploited in Bech and Andersson (1996b) 
to quantify the structure of the Reynolds stress in the central region of Couette flows. For 
Couette flows Eq. (18) gives vm ~ ^Vx/4:0. The simulation just mentioned has a Reynolds 
number^^ of 5200, i.e. well above the threshold of transition to turbulence. The behavior 
of the Reynolds stress as a function of the distance to the wall is represented on Fig. la of 
Bech and Andersson (1996b), and, after accounting for the particular normalization adopted 
in their graph^^, one flnds vm — for this simulation, which is nearly identical to the 

estimate deduced from Eq. (18). Even if one takes into account the fact that the value of 

which is relevant for the bulk of the flow is smaller than the one adopted in Eq. (18) by 
a factor ~ 4, the two estimates of vm still agree within a factor of ~ 3. Another estimate 



^^Note that our definition of the Reynolds number differs from the one adopted in these papers by a factor 
of 4. 

^^A property of Couette flows is that the total mean shear stress r = p{jyd{vx) / dy — (Sv^Svy)) is constant 
in the shearwise direction (this follows from the stationarity of the mean flow). Away from the wall t ~ 
—p{5vxSvy) whereas close to the wall r ~ pvd{vx) / dy; consequently the Reynolds stress is usually normalized 
to r/p, velocities to Vr = (r/p)^/^, and on has Vm — fr in the bulk of the flow. The value of Vr for this 
simulation can be obtained in the following way. Note first that Fig. (la) of their paper also displays 
Re^^d{vx) /dy where, for their simulation, Rct = Vrh/u = 82, and where Vx is normalized to Vt and y to the 
walls half-distance h; this quantity is equal to 1 in the immediate vicinity of the wall. On the other hand, 
the value of the velocity gradient near the wall in units of 2Uw/h = Avx/h can be deduced from Fig. 4 of 
Bech et al. (1995) which relates to the same simulation. The comparison of these two measures of the same 
quantity yields the required value of Vr/2Uyj. 
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of the same quantity for non-rotating free shear flows is obtained from the representation^^ 
of (Sv^Svy) in Fig. 14 (at Q = 0) of Bidokhti and Tritton (1992), and gives vm ~ Afj-ZlO, 
which differs from the order of magnitude estimate quoted above by a factor of ~ 4. Some 
indication of the value of Im can also be extracted from Fig. 2 of Bcch et al. (1995), which 
shows the power spectra in the shearwise and spanwisc directions. However, the box size in 
these directions is large in terms of /i, and the resolution of the simulation does not allow 
the authors to really reach the inertial part of the turbulent spectrum. This is particularly 
noticeable for spectra in the middle of the flow (displayed in the y = 82 quadrant of this 
figure), which are nearly fiat down to kxh ~ 10, and drop precipitously for larger values 
of k because of numerical dissipation, as the limit resolution of the simulation is reached. 
The kz spectra behave sensibly better, more probably because the box is 2.5 times smaller 
in this direction, and show some indications that an inertial spectrum tries to develop for 
8 ^ kzh < 30. Because h = Ay/2, this suggests that Im in this simulation is at most 
within a factor of ~ 3 of the order of magnitude estimate deduced from Eq. (17). Note in 
passing that, for Couette flows, the inertial spectrum docs not need to be resolved in order 
for turbulence to be observed in numerical simulations; this is related to the existence of a 
large scale nonlinear mechanism which sustains turbulence, as mentioned in section 2.1, and 
which is most likely at the origin of the more or less fiat part of the spectra at large scales. 

REFERENCES 

Andereck, CD., Liu, S.S., and Swinney, H.L. 1986, J. Fluid Mech., 164, 155 

Balbus, S.A., and Hawley, J.F. 1991, ApJ, 376, 214 

Balbus, S.A., Hawley, J.F., and Stone, J.M. 1996, ApJ, 467, 76 

Bech, K.H., Tillmark, N., Alfredsson, P.H., and Andersson, H.I. 1995, J. Fluid Mech., 286, 
291 

Bech, K.H., and Andersson, H.I. 1996, J. Fluid Mech., 317, 195 
Bech, K.H., and Andersson, H.I. 1996, Fluid Dyn. Res., 18,65 
Bech, K.H., and Andersson, H.I. 1997, J. Fluid Mech., 347, 289 
Bidokhti, A.A., and Tritton, D.J. 1992, J. Fluid Mech., 241, 469 

^•^Note that vm characterizes turbulent transport, and that {SvxSvy) is usually smaller than the magnitude 
of velocity fluctuations by a factor of a few. 



Brandenburg, A., Nordlung, A., Stein, R.F., and Torkelsson, U. 1995, ApJ, 446, 741 
Casse, F., and Ferreira, J. 2000, A&A, 353, 1115 
Chandrasekhar, S. 1960, Proc. Nat. Acad. Sci., 46, 53 
Dauchot, O., and Daviaud, F. 1994, Phys. Fluids, 7, 335 
DubruUe, B. 1993, Icarus, 106, 59 

DubruUe, B., and Zahn, J.-P. 1991, J. Fluid Mech., 231, 561 
Ferreira, J. 1997, A&A, 319, 340 

Fleming, T.P., Stone, J.M., and Hawley, J.F. 2000, ApJ, 530, 464 
Johnson, J.A. 1963, J. Fluid Mech., 17, 337 

Hamilton, J.H., Kim, J., and Waleffe, F. 1995, J. Fluid Mech., 287, 317 
Hawley, J.F., Gammie, C.F., and Balbus, S.A. 1995, ApJ, 440, 742 
Hawley, J.F., Balbus, S.A., and Winters 1999, ApJ, 518, 394 
Hawley, J.F., Balbus, S.A., and Stone, J.M. 2001, ApJ, 554, L49 
Hur, J.-M., Richard, D., and Zahn, J.-P. 2001, A&A, 367, 1087 
Jimenez, J., and Moin P. 1991, J. Fluid Mech., 225, 213 

Komminaho, J., Lundbladh, A., and Johansson, A.V. 1996, J. Fluid Mech., 320, 259 
Landau, L.D., and Lifshitz, E.M. 1987, Fluid Mechanics, 2nd edition (Pergamon Press) 
Lesieur, M. 1987, Turbulence in Fluids (Kluwer Academic Press) 
Lerner, J., and Knobloch, E. 1988, J. Fluid Mech., 189, 117 
Marcus, P.S. 1984, J. Fluid Mech., 146, 45 
Marcus, P.S. 1984, J. Fluid Mech., 146, 65 
Porter, D.H., and Woodward, P.R., ApJS, 127, 159 
Prandtl, Z.A. 1925, Zs Angew. Math. Mech., 5, 136 



-28- 



Richard, D. 2001. Instabilites hydrodynamiques dans les ecoulements en rotation 
differentielle. PhD thesis. Universite de Paris VII. 

Richard, D., and Zahn, J.-P. 1999, A&A, 347, 734 

Sano, T., Miyama, S.M., Umebayashi, T., Nakano, T. 2000, ApJ, 543, 486 
Taylor, G.I. 1936, Proc. Roy. Soc. London, A 223, 289 

Tennekes, H, and Lumley J.L 1972, A first course in turbulence (Cambridge: MIT Press) 
Terquem, C, 2001, astroph/0107408 

Tillmark, N., and Alfredsson, P.H. 1992, J. Fluid Mech., 235, 89 
Tritton, D.J. 1992, J. Fluid Mech., 241, 503 
Waleffc, F. 1997, Phys. Fluids, 9, 883 
Wendt, F. 1933, Ing. Arch., 4, 577 

Zeldovich, Y.B. 1981, Proc. Roy. Soc. London A, 374, 299 



This preprint was prepared with the AAS IM^jX macros v5.0. 



- 29 - 



u 

»- 




Fig. 1. — Sketch of the configuration of Couette flows. The flow is bounded by two counter- 
moving walls, and boundary layers develop in the turbulent regime, as shown by the mean velocity 
profile. By putting the experimental setup on a rotating platform, one obtains the so-called rotating 
Couette flow. 
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Fig. 2. — Sketch of the configuration of (rotating) free shear layers. Two layers of fluid of different 
velocities, initially horizontally separated, come in contact at the end of a dividing plate, and a 
turbulent shear layer develops and widens downstream. 



- 31 - 




Fig. 3. — Idealized behavior of the minimal Reynolds number of fully turbulent Couette- Taylor 
flows, function of the relative gap width. 



